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04 ■ Abstract 

As an image of the many-to-one map of loop-erasing operation £ of random walks, a self-avoiding 
walk (SAW) is obtained. The loop-erased random walk (LERW) model is the statistical ensemble of 
SAWs such that the weight of each SAW £ is given by the total weight of all random walks n which 

Si ' 

are inverse images of £, {it : C(tt) = £}. We regard the Brownian paths as the continuum limits of 
random walks and consider the statistical ensemble of loop-erased Brownian paths (LEBPs) as the 

a . 

continuum limits of the LERW model. Following the theory of Fomin on nonintersecting LERWs, 
we introduce a nonintersecting system of A-tuples of LEBPs in a domain D in the complex plane, 
where the total weight of nonintersecting LEBPs is given by Fomin's determinant of an N x N 



OS 



CO 

> 

matrix whose entries are boundary Poisson kernels in D. We set a sequence of chambers in a planar 
domain and observe the first passage points at which N Brownian paths (71, . . . , 7jv) first enter 
each chamber, under the condition that the loop-erased parts (£(71), . . . , £(7^)) make a system 
of nonintersecting LEBPs in the domain in the sense of Fomin. We prove that the correlation 
functions of first passage points of the Brownian paths of the present system are generally given 
by determinants specified by a continuous function called the correlation kernel. The correlation 
kernel is of Eynard-Mehta type, which has appeared in two-matrix models and time-dependent 
matrix models studied in random matrix theory. Conformal covariance of correlation functions is 
demonstrated. 
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I. INTRODUCTION 



The vicious walker model introduced by Fisher [II] is a one-dimensional system of simple 
symmetric random walks conditioned so that any pairs of trajectories of walkers in the 1+1 
dimensional spatio-temporal lattice are nonintersecting. The transition probability of iV vi- 
cious walkers can be described by using a determinant of an N x N matrix, whose entries are 
transition probabilities of a single random walker with different initial and final positions. 
This determinantal expression for nonintersecting paths is called the Karlin-McGregor for- 
mula in probability theory j^J and the Lindstrom-Gessel-Viennot formula in enumerative 
combinatorics s[ Preserving the determinantal expression for transition probability den- 
sities, a continuum limit (the diffusion scaling limit) of vicious walkers can be taken and 
we have the system of one-dimensional Brownian motions conditioned never to collide with 
each other 5KZ|. The important fact is that the obtained interacting particle systems de- 
fined in the continuous spatio-temporal plane, which can be called the noncolliding Brownian 
motion is identified with Dyson's Brownian motion model with {3 = 2 jjj, which was 
originally introduced as a stochastic process of eigenvalues of an Hermitian-matrix valued 



Brownian motion in the random matrix theory [10|, 111] . The notion of correspondence be- 
tween nonequilibrium particle systems and random matrix theories is very useful 12] and 
spatio-temporal correlation functions of noncolliding diffusion processes have been deter- 
mined explicitly not only for the systems with finite numbers of particles but also for the 



systems with infinite numbers of particles 



ay 



16]. 



In the present paper, we study a system of continuous paths not in the 1+1 dimensional 
spatio-temporal plane but in the two dimensional plane (i.e. the complex plane C = {z = 
x + iy} with i = y/— 1), which will be called the nonintersecting system of loop-erased 
Brownian paths (LEBPs) [l7]. A version of nonintersection condition is imposed between the 
)at hs ( see Eq. (J7J) below) , and then the total weight of LEBPs is given by Fomin's determinant 



18, 



191 ] instead of the determinant of Karlin-McGregor (and of Lindstrom-Gessel-Viennot). 



There the entries of matrix whose determinant is considered are the normal derivatives at 
boundary points of domain of the Green's functions of the two-dimensional Poisson equation 
(the Poisson kernels and the boundary Poisson kernels) instead of the transition probability 



densities 



In Section II, we define the LERW model and briefly review Fomin's theory of noninter- 
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secting LERWs. 

As the continuum limit of LERW model, the statistical ensemble of LEBPs is introduced 
in Section III. A. There the Green's function, the Poisson kernel, and the boundary Poisson 
kernel are defined for the Brownian motion in a domain in the two-dimensional plane or the 
complex plane C. Then for an L x ir rectangular domain Rl in C, Fomin's determinant of 
the boundary Poisson kernels and of the Poisson kernels for iV-tuples of Brownian paths are 
studied, and nonintersecting system of LEBPs are constructed in Section III.B (see Figj3]). 
In Section III.C, we set two rectangular domains on C adjacent to each other at a vertical 
line Re z = x, in which N Brownian paths are running from the left rectangular domain 
to the right one through the line Kez = x (see FigS]). We impose the condition that the 
loop-erased parts of Brownian paths are nonintersecting in the sense of Fomin as expressed 
by Eq.([7]). Under this condition, the probability density function of the first passage points 
on the line Kez = x, at which the N Brownian paths enter the right rectangular domain 
from the left domain, are given as Eqs.( |36|) for L < oo and ( |39l) for L — > oo, respectively 
In Section III.D, we consider a sequence of M + 1 rectangular domains on C, M e N = 
{1,2,...}, where the m-th domain and the (m + l)-th domain are adjacent at the line 
Kez = x m , 1 < m < M. iV-tuples of Brownian paths are running from the left to the right 
(see Figj5]) under the condition that their loop-erased parts make a nonintersecting system 
of LEBPs. The probability density function of joint distributions of first passage points at 
M lines Rez = x m , 1 < m < M, of the Brownian paths are determined as Eqs. (l43l for 
L < oo and (jSj) for L — > oo, respectively. 

In Section IV, a special initial condition is assumed when the iV-tuples of Brownian paths 
start from the left boundary of the leftmost domain. In this special case, we can explicitly 
obtain all multiple correlation functions of first passage points on the lines Kez = x m , 1 < 
m < M for any M G N, in which they are given by determinants (Theorem 1). The 
correlation kernel, which completely specifies the determinants, are given by Eqs.( |54l) and 
(1551) . The statistical ensemble of points whose correlation functions are generally expressed 
by determinants with a correlation kernel is called a determinantal point process or a Fermion 



point process in probability theory |2Q|-|22j|. It should be noted that the present correlation 
kernel is asymmetric K 7 ^ 2 (x,8; x' ,6') ^ K^ 2 (x', 8'; x, 9) for x ^ x' as shown by Eqs. (l54|) 
and fl55l) . This asymmetric correlation kernel is of Eynard-Mehta type j^, which has been 
studied in two-matrix models and time-dependent matrix models in random matrix theory 
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[5. 



sua mi. 

Since the Brownian motions and their loop-erased parts on C are conformally invari- 
ant [17], our correlation kernel is conformally covariant. In Section V, the determinantal 
correlation functions in the half-infinite-strip domain R = UmL^ooRi = {z G C : Rez > 

0. < Imz < tt} given by Theorem 1 is mapped to the domain Q = {z = re ld G C : r > 

1, < 9 < 7t} (Corollary 2). There numerical plots of the density function and the two-point 
corrrelation functions in the domain Q are shown by figures. 

Concluding remarks are given in Section VI. Appendix A is prepared to derive the for- 
mulas of the Poisson kernel and the boundary Poisson kernel used in the text. 



II. LOOP-ERASED RANDOM WALKS AND FOMIN'S DETERMINANT 

We consider an undirected planar lattice consisting of a set of vertices (sites) V = {vj} 
and a set of edges (bonds) E = {ej}. Together with a set of the weight functions of the 
edges W = {w(e)} e€ E, a network T = (V, E, W) is defined. 

For a, b G V, let tt be a walk given by 

7T : a — Vq — Vt> i -4-i>2 ■ • ■ v m — b (1) 

where the length of walk is \ir\ = m G N and, for each < j < m — 1, Vj and Vj+i are 
nearest-neighboring vertices in V and ej G E is the edge connecting these two vertices. We 
will shorten (pQ) to n : a — » 6, or a A b. The weight of n is given by w(n) = YYjLi w { e j)- For 
any two vertices of a, b G V, the Green's function of walks {n : a — ► b} is defined by 

7r:a— ¥b,\ir\=m 

The matrix VT^ = (W(a, 6)) a ,6ev is called the watt matrix of the network T. 

The loop-erased part of tt, denoted by >C(7r), is defined recursively as follows. If 7r does 
not have self-intersections, that is, all vertices Vj,0 < j < m are distinct, then C(n) = it. 
Otherwise, set C(tt) = C(tt'), where tt' is obtained by removing the first loop it makes. In 
other words, if (k, £), k < £ is the smallest pair in the index set {j '• : < j < m} of {v jj^Lo * n 
the sequence ([T]) such that = Vg, then the subsequence ^ v k+1 A- 2 • • • is removed 
from 7T to obtain tt'. 
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The loop-erasing operator £ maps arbitrary walks to 'walks without self-intersections', 
which are usually called self-avoiding walks (SAWs). Note that the map is many-to-one; if 
C is a SAW obtained by applying C, the set of inverse images {n : C{ji) = (} has more than 
one element in general. For each SAW (, the weight w(() is given by 

w(o= Yl ( 3 ) 

7r:£(7r)=C 

We consider the statistical ensemble of SAWs with the weight (EJ) and call it loop-erased 
random walks (LERWs) |24j. Note that the LERW model is different from the SAW 
model, in the sense that, though the configuration space of walks are the same, the weight 
of each walk is different from each other. In the SAW model we consider a statistical 
ensemble of SAWs with the weight w(Q = e _/3 '^', where e 13 is the SAW connective constant, 
while the weight of SAW in the LERW model is given by the sum of weights of all walks, 
which are the inverse images of the projection C as shown by (j3J). 

Assume that A = (ai, a 2 , . . . , a/v) C V and B = (b\, b 2) . . . , 6/v) C V are chosen so that 
any walk from a,j to bk intersects any walk from aj',j' > j, to bk>,k' < k. The weight of 
iV-tuples of independent walks a\ 61, . . . , is given by the product of iV weights 

Y[f=i w ( 7l e)- Then we consider iV-tuples of walks (711, 7r 2 , . . . , 7rjv) conditioned so that, for 
any 1 < j < k < N, the walk 7Tk has no common vertices with the loop-erased part of 7Tj] 

Cfa) n 7r fe = 0, 1 < j < k < N. (4) 

See FigJTJ By definition, C(iTk) is a part of 7?^, and thus nonintersection of any pair of 
loop-erased parts is concluded from (TJJ; 

£(7Tj) n £(7r fc ) = 0, 1 < j < k < N. (5) 

Fomin proved that total weight of iV-tuples of walks satisfying such a version of nonin- 
tersection condition is given by the minor of walk matrix, det(W / A,B) = det a< =A,beB(W(a, b)) 
1^1 . This minor is called Fomin's determinant in the present paper and Fomin's formula is 



expressed by the equality [17H19] 



N 

det(W AB )= H W M- ( 6 ) 

C(ivj)mt h =<t,j<k i=i 

It should be noted that the RHS of (jSJ) seems to depend nontrivially on the ordering of the 
sets A and B, but it is indeed invariant up to a sign under the change of ordering, since it is 




FIG. 1: The situation C(iTj) D tt^ = 0,j = 1,2 is illustrated in a planar domain D, where A = 
(01,02,03) and B = (61,62,63) are all boundary points of dD. In this figure, C{tti) and £(^2) 
denoted by solid curves are the loop-erased parts of the walks n\ : a\ — )■ 61 and TT2 : 02 — )• 62, 
respectively, where loops described by broken curves are erased. The third walk 7r3 : 03 — > 63 can 
be self-intersecting, but it does not intersect with C(tti) nor £(^2). As a matter of course, £(^3) 
is a part of ^3, and thus £(7Tj) H £(^3) = 0, j = 1, 2. 

equal to the LHS of ([6]), which is antisymmetric in exchanging any pair of rows or columns 
18]. 

III. NONINTERSECTING SYSTEMS OF LOOP-ERASED BROWNIAN PATHS 
A. Fomin's determinants for Brownian paths 

Random walks on a lattice can be regarded as discrete approximations of a Brownian 
path. For example, for a two-dimensional Brownian path 7 starting from a point a and 
terminating at a point 6 in a domain D in the plane R 2 , we can consider a random walk on 
a planar lattice T embedded in R 2 to approximate 7. Assume that the lattice spacing of the 
planar lattice is 1/n, n G N, choose two vertices and such that they are nearest 
vertices in to a and b, respectively, and denote the random walk on from to 
M n ) by n( n \ The continuum limit can be taken by getting n — > 00 and ir^ will converge to 
a Brownian path 7 running from a to 6. When the points a and b are on the boundary of a 
domain and especially when the boundary is not smooth, we need careful consideration 
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for continuum limit even for a single path. See 25] for rigorous argument. 

Then, if we can perform the same limiting procedure n — > oo for Fomin's determinants 
of iV-tuples of random walks {ir[ n \ . . . , vrjv^} on planar lattices r^, the limit will give the 
total weight of the iV-tuples of two-dimensional Brownian paths (71, . . . , 7jv), which satisfy 
the following condition, 

£(7i) n lk = 0, 1 < 3 < k < N, (7) 

where £(7?) is considered to be the n — > 00 limit of the sequence of loop-erased parts 
neN of discrete approximations {7rj n ^} ne N of 7,, 1 < j < N. 
Recently Kozdron and Lawler [it| gave the mathematical justification of the above men- 
tioned continuum- limit-procedure for simply connected planar domains D in the complex 
plane C, where the initial and the final points A = {a,j} and B = {bj} of paths {7,} can 
be put on the boundaries of the domains 3D. We should note that the characteristics of 
Brownian motion look more similar to those of a surface than those of a curve (see, for 
example, Chapter 1 of 26|] ) . It implies that the Brownian path has loops on every scale and 
then the loop-erasing procedure defined for random walks in Sect. II does not make sense for 
Brownian motion, since we can not decide which loop is the first one. Kozdron and Lawler 
proved explicitly, however, that the continuum limit of Fomin's determinant of the Green's 
unctions of random walks converges to that of the Green's functions of Brownian motions 



171 ] . This will enable us to discuss nonintersecting systems of LEBPs without dealing with 



individual LEBP. See a remark put at the end of the present paper and 13, ]27J, 128] for more 
details. 

One of the advantages of taking continuum limit is the fact that the Green's function 
and its normal derivatives at boundary points (the Poisson kernel and the boundary Poisson 
cernel) are obtained by solving the Laplace equation with appropriate boundary conditions 



171 ] . As a realization of the two-dimensional Brownian motion, we consider a complex 



Brownian motion B t = B^ + iBf \ where B^ and B^ are independent one-dimensional 
standard Brownian motions satisfying (dB^) 2 = {dBf ) 2 = dt and dB^dBf = 0, t > 0. 
For a domain D C C, let pD{t,a,b) be the transition probability density of the complex 
Brownian motion from a G D to b G D with duration t > with the absorbing boundary 
condition at 3D. The Green's function for this Brownian motion is defined by 

POO 

Gn(a,b) = / p D (t 1 a 1 b)dt. (8) 
Jo 



(Note that the summation with respect to the length of walk m = \tt\ in (121) for the Green's 
function of random walks is here replaced by the integral with respect to duration of time 
of Brownian motion.) It is also the Green's function for the Poisson equation A a Gc(a, b) = 
S(a—b), where A a is the Laplacian with respect to the variable a, with the Dirichlet boundary 
condition Gr>(a,b) = on a G 3D. The complex Brownian motion is conformally invariant 
in the sense that the Green's function in a domain Gd(ci, b),a,b G D has the property 

G D (a,b) = G D ,(f(a)J(b)) (9) 

for any conformal transformation / : D — > D' . Then we will select a suitable domain D in 
C such that the Poisson equation can be analytically solved and explicitly determine the 
Green's function Go, and then, though the equality (J9]), we can obtain Gd' for other domain 
D' by an appropriate conformal transformation. 

For a G D and b G 3D, the Poisson kernel Hjj(a, b) is defined by 

H D (a, b) = \\iva -G D (a, b + en b ), (10) 

where denotes the inward unit normal vector at b G 3D. By this definition, we see 
that Ho(a,b) solves the Laplace equation A a Ho(a,b) = for a G D, that is, Ha(a,b) is a 
harmonic function of a G D, and satisfies the boundary condition 

lim H D (a,b) = 6 (a - b) for a G dD. (11) 

a— >a:a&D 

Moreover, we define the boundary Poisson kernel Hgi)(a, b) for a, b G dD by 

H 9D (a, b) = lim -H D (a + en a , b). (12) 

From the conformal invariance ([9]) of the Green's function Gd, the following conformal 
covariance properties are derived; if / : D — > D' is a conformal transformation, 

HD(a,b) = \f(b)\H D '(f(a)J{b)), a eD.be 3D, (13) 
H dD (a,b) = \f'(a)\\f(b)\H B D.(f(a)J(b)), a,be3D. (14) 



See Section 2.6 of [17| and Chapters 2 and 5 of [29] for further information about Poisson 
kernels, boundary Poisson kernels and their conformal covariance. We will study Fomin's 
determinants of the form deti<^fc<jv and deti<j^<N[HdD{.0'j, b k )] in the following. 



S 



B. System in a rectangular domain and the crossing exponent 



In order to give explicit expressions for Poisson kernel and boundary Poisson kernel, now 
we fix the domain D as the following rectangular domain 

R L = {2 = x + iy G C : < x < L,0 < y < 7r|, L > 0. (15) 

As shown in Appendix A, for < x < L, < 9, p, ip < it, the Poisson kernel connecting an 
inner point x + \9 G D and a point L + ip at the right boundary of Rl, dR^ = {L + iy : < 
y < 7r}, is given by 

_ _ > .„ , . . 2 *^ sinh(nx) sin(ra#) sin(np) . . 

F K£ 3 + iS, L + ip = - 1 . rr , 16 

7r ^— ' sinn nL 

n=l v ' 

and the boundary Poisson kernel connecting a boundary point i<^ on the left boundary of 
Rl, dR\ = {iy : < y < n}, and a boundary point L + ip on the right boundary dR^ is 
given by 

.. _ . . 2 ^ nsin(no9) sin(np) , s 

H aRL (i<p, L + ip) = -J2 . I n ■ 17 

ii ^-^ sinn nL 

n=l 

The definitions (jH}, ffTUl) and f|T2l imply that H Rl (x + i9, L + ip)dp gives the total weight 
of the paths of complex Brownian motion starting from the inner point x + 16 £ Rl, which 
make first exit from the domain Rl at a point in the interval [L + ip, L + i(p + dp)] on the 
right boundary dR^, and that Hq Rl (iip,L + ip)dipdp gives the total weight of the paths of 
complex Brownian motion, which enter Rl at a point in [i<p, i(<p + dip)] on dR\, and make 
first exit from Rl at a point in [L + ip, L + i(p + dp)] on dR^. (See Figj2j) 

By the conservation of probability, the following equalities should be satisfied; for < 
x < x' < L, 

H RL (x + i6,L + ip) = / d6' H R ,(x + i6,x' + i6')H RL {x' + i6',L + ip), (18) 

Jo 

H dRL (iip,L + ip) = / d0" H aRx (i<p,x + W")H Rli (x + iB",L + ip), (19) 
Jo 

< ip, p < 71. As a matter of course, the validity of them can be directly confirmed by apply- 
ing the orthogonality relation of the sine functions, f Q d9 sm(n8) sm(m8) = (tt /2)S nm , n,m G 
N, to the expressions (IT5|) and (IT7|) . We note that in Eq. (IT8l) the point x' + 16' is regarded 
as the first passage point on the line Re z = x' in Rl of the Brownian path running from 
x + 16 G Rl to L + ip G and similarly that in Eq.(fl9~]) a; + i6*" is regarded as the first 




FIG. 2: The broken curve denotes a path of complex Brownian motion starting from an inner point 
x+i6 E Rl, which makes first exit from Rl at L+ip E dR^. The solid curve does a path of complex 
Brownian motion, which enters Rl at 'up E R\ and makes first exit from Rl at L + ip E dR^. The 
former path contributes to Hr l (x + W, L + ip) and the latter does to Hqr l (i<p, L + ip), respectively. 
The point x' + i9' is the first passage point on the line Re z = x' in Rl of the former path, and the 
point x + iO" is the first passage point on the line Re z = x in Rl of the latter path. 

passage point on the line Kez = x in Rl of the Brownian path running from 'up E dR\ to 
L + ip E See Figf5J 

For N e N, let = {0 = (9 U 9 2 , ■ ■ ■ , 6 N ) : < 9 X < 9 2 < • • • < 9 N < vr}. Then, for 
cp = (ifi, . . . , ip N ) E Wjy-, p = (pi, . . . , p^) E W]y, consider Fomin's determinant 

(20) 

If we write the path of j-th Brownian motion, 1 < j < N, which enters Rl at 'ufj E dR\ 
and makes first exit from Rl at L + ipj E dR^ as jj, then ( |20l) gives the total weight of the 
TV-tuples of paths (71, . . . , 7^) satisfying the nonintersection condition with the loop-erased 
parts ([7]). We note again that it is a sufficient condition for 

£(7,) n C( lk ) = 0, l<j<k<N. (21) 

See FigEl 



HdR L (i<Pj,L + ip fc ) 
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FIG. 3: The condition ([7]) is illustrated for k = 3, N = 3. The Brownian path 73 : itps — > L + ip3 in 
Rl can be self-intersecting, but it does not intersect with £(71) : i</?i — >• L + ipi nor £(72) : i<^2 — > 
L + ip2- By definition £(73) is a part of 73, and thus £(71) H £(73) = and £(72) H £(73) = 0. 

By multilinearlity of the determinant, we find that Eq. f[2"Uj) with Eq.( lT7|) is written as 



iV 



e n 



L sinh(n,L) 

n=(ni,...,n JV )GN JV 3=1 V J 7 

>V! 2 1< det jv [sin(n CT0 - ) ^)sin(n CT0) p fc ; 



N 



det sin(jV&; 

7T / l<j,k<N 1 



det 

J K£,m<N 



sm(£p m )] x^2a x sx((fi)s x (p), (22) 



where Sjv is a set of all permutations {a} of {1, 2, . . . , N}, X = (Ai, . . . , A./v) with Aj 
njv-j+i - (iV - j + 1), 1 < j < iV, and 

+ iV - j + 1 



AT 

n 



L S inh((A i + ^-j' + l) J L)' 

sin((A fe + iV- fc + l)^ 



det 

l<j,k<N 



det 

l<j,k<N 



sin((iV- k + l)ipj) 



(23) 



(This is a modified version of 'Schur function expansion' used in [8|, [l2j.) We note that 

l(e-i)/2] 



sm(£9) = sin 9 



2 £ - 1 ('cosi 



- + E i-^T 1 ) 



(2cos#) 



!-2s-l 



(24) 



where, for r 6 R, [r] denotes the greatest integer not greater than r. Then for 6 



det 

Kt,m<N 



sin 



det 

l<e,m<N 



2 e 1 sin 9 m \ (cos 9, 



0((cos9 r 



2 N(N-l)/2j lN{e) 



(25) 
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with the function 



N 



sin 6j J f (cos 6e - cos 9 k ). 

j=l \<k<l<N 



(26) 



Since 



N 



2 N H(\ j + N-j + l)e- L ^U^+^+V as l^oo, 



and in particular for = (0, 0, . . . , 0) 

a^2 N N\e~ LN ^ N+1 ^ 2 as L ->• oo, 
we can conclude from the expansion (1221) that 

/» (L, p|<p) ~ 2 N ^^ N N\ e- N ^ L l 2 h N ^)h N {p) as L ^ oo. (27) 
Since the simple iV-product of the boundary Poisson kernels Ylj=i HdR L (i<Pj, L + ipj) ~ 



nf =1 2/{7rsinhL} ~ (4/ir) N e- NL as L — >• oo, the ratio behaves 



/£(£,p|v) 



A RL (ip,p) = 



c N (cp,p)e 



-ip N L 



as L — )■ oo 



with 



^ N =-N(N-1), 



(28) 



(29) 



where cjv(</?, p) = 2 N<yN ^ Nlhw ((f) fiN (p) • The exponent (|29|) is called the crossing exponent 



C. Probability density function for first passage points of Brownian paths under- 
lying in nonintersecting LEBPs 

We consider the integral of Fomin's determinant (|20l) over all possible ordered sets of 
exits p G W^r at the right boundary dR^, 



dpf N (L, P \v) 



VV JV 



dp det 

W - l<j,k<N 

vv AT 



H d R L (i(pj,L + ip k ) 



\ dp det / d0# 9i k(i^j,x + i0)#tfja; + i6^ + ip fc ) 
y W x i<j,fc<iv y 



(30) 



where dp = Ylj=i dpj, < x < L, and Eq. (fT9|) has been used in the last equality. Applying 



the Heine identity 



dz det 

l<j,k<N 



4>j(z k ) 



det 

K£,m<N 



ipt{zm)] = det 



dz(pj(z)tlj k (z) 



(31) 
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which is valid for square integrable functions <pj,ipj, 1 < j < N, it is written as 



dp 



d0 det 

l<j,k<N 



H aRx (i(pj,x + \9 k 



det 

KLm<N 



H Rl (x + i9i,L + ip n 



dO fx(x,0\ip) / <ip det H Rl (x + i9g, L + ip 



l<e,m<N 



(32) 



Then, if we introduce the integral 



Af N (x,L,6) = dpf N {x,6;L,p), 



< x < L, 



(33) 



of Fomin's determinant for the Poisson kernels 



det 

l<j,k<N 



f N (x,0;L,p 
6, p e Wjy, and divide the both sides of 

I / d0f°(x,0\cp) 



H Rl (x + i9j,L + ipk) 



(34) 



by Af^(L, (p), we obtain the equality 



(35) 



Then, given ip e W]y-, if we put 



Af N (x,L,6) 



< x < L, 



(36) 



Nn(L, cp) ' 

for G W^r, it can be regarded as the probability density function. As illustrated in FigJH 
here we consider the rectangular domain (fl5l) as a union of the two rectangular domains 
R x = {z e C : < Rez < x,0 < lmz < ir} and R{ x ,l} = {z E C : x < Re z < L,0 < 
Imz < 7r}, which are adjacent to each other at a line Rez = x. We have considered iV-tuples 
of Brownian paths (71, ... , 7^) starting from (iifi, . . . , i<Pn), all of which run inside of the 
domain Rl = R x U R{ x ,l} until arriving at the right boundary dR^, under the condition 
that (£(71), . . . , £(7jv)) makes a system of nonintersecting LEBPs in the sense of Fomin (JTj). 
In an ensemble of such iV-tuples of Brownian paths, Eq. fl36|) gives the probability density 
of the event such that the first passage points of (7i,...,7jy) on the line Rez = x are 
(x + . . . , x + i9n). That is, for 1 < j < N, 7^ makes a first exit from the left rectangular 
domain R x and enters the right rectangular domain R{ x ,l} at x + i9j. The path jj, which 
made a first exit from R x at x + i9j, can reenter R x at different point x + \9' on the line 
Rez = x. As shown by the path 72 starting from itp2 in FigJU the path jj can make a loop 
passing the line Rez = x, and if the first passage point x+i9j is included in such a loop, that 
point can not be included in C(jj). In other words, Eq. fl36|) gives the probability density for 
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x L 



FIG. 4: Three Brownian paths (71,72,73) in Rl = U R{ X ,L} starting from (i(pi,'vp2,'vp3) and 
arriving at dR^ whose loop-erased parts (£ (71), £(72), £(73)) are nonintersecting in the sense of 
Fomin. The first passage points on the line Re z = x are {x + i9i,x + 162, x + 1O3). Since x + 162 is 
in a loop of 72 as shown by a broken curve, it is not included in £(72). The probability density of 
points (x + Wi, x -\- i02, x + 16*3) is given by Eq. (f36j) . 

N nonintersecting loop-erased Brownian paths (£(71), . . . , £(7^)) starting from the points 
(i(fi, ■ ■ ■ , i<^jv) and satisfying the condition ([7j) in R L , such that the underlying Brownian 
paths (71, . . . , 7jv) are realized so that they first arrive at the vertical line Re 2; = x < L at 
(x + i0 1 , . . . , x + i9 N ). 

Now we consider the system in the limit L — > 00. By applying the asymptotics (1271) . we 
have 





2 N ( N +^n- N N\e- N ^ N+1 ^ 2 h N (cp) x 




dp h N (p) 



(37) 



as L 00. Similarly, we can obtain the following asymptotics of ( |33l) 



f 2 °° 

Af N (x,L,0)= / dp det -> 

^ vv iv L n=l 




^2 N ^ N+1 ^~ N e- N ^ N+1 W 2 h N (e)l[smhUx) / dp^(p), 




(38) 
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as L — > oo. Then, given cp G W]y, for G W 



7T 

N ■ 



where 



p N (x,0\<p) = lim p^(x,0|c^>) 



CW(s)/&(x,0|¥O 



A/ 



, < a; < oo, 



(39) 



(40) 



D. Joint distribution of first passage points in a sequence of chambers 

Let M G N and < x\ < x 2 < ■ ■ ■ < xm < L < oo. Here we consider M vertical 
lines in R N at Rez = x m , 1 < m < M. For (x m ,x m+1 ), (m) = (^ m) ,...,^ } ) G W n N , 
0(m+i) = ^(m+i)^ ^ (m+i)^ ^ W^, 1 < m < M — 1, we define 

^(z m ,0^;x m+1 ,0( m+1 >) = A/v(awi, L, ^ 1) ) /jy(a - m> M gfrn+i)), (41) 

A/^(x m ,L,0 (m) ) 

Then by definition (|36|) . we have 



d0« ^(xi, (2) ) 

I? 

w- Nn(L,<P) Af N ( Xl ,L,0 (1 >) 



fe^P / do^f N ( Xl ,e^)f N ( Xl ,e^; X2 ,e^). 

By the Heine identity (l3Tj) and the equality (TT9~j) . 

^^(xi,^ 1 ^)/^^,^;^,^) 



d0 (1) det 



k / 



det 

J l<e,m<N 



H R (x 1+ i9f\x 2 + ie$) 



det 

l<j,k<N 



det 

l<j,k<N 



and thus we obtain the equality 



l<j,k<N 

dB® H dRxi (iifj, x x + i6^)H Rx2 ( Xl + i#«, x 2 + i0 



(2b 



15 



Obviously this equality can be generalized as 

p L N (x m+ i,0 {m+1) \v) = [ rf0 (m) ^(a; m ,0 (m V)^(^,0 (m) ;x m+1 ,0 (m+1) ) (42) 

Jw- N 

for 1 < m < M - 1. 

Then, given ip G W^, if we introduce a function of (m) = (6 { ™\ . . . , 0#°) G Wfr, 1 < 
m < M by 



M-l 



Pv 



On,0 (1) |¥>) II 4(^,0 (m) ;^ m+ i^ (m+1) ) 



m=l 



Kf I T q{M)\ M—l 

= \fd(j \ ] f d N^O {1) \v) II fN(x m ,6^; Xm+1 ,e^), (43) 
it can be regarded as the probability density function, since it is well normalized as 
IT / dOWp L N (x 1 ,OM;...;x M ,OM\<p)= [ d0™ p L N (x M ,e^\cp) = 1 

m=l ^ vv jv ^ vv jv 

by ( 142)) and fl35|) . As shown by Figj5j here we think that Rl has M + 1 'chambers' 

#{* ra -i,* m }, 1 < m < Rl = Um=i R {x m -i,x m } with x o = °> = L, where i^.^} 

and -R{ Xm ,a: m+ i} are adjacent to each other at the line Rez = x m , 1 < m < M. Under the 
initial condition that jj starts from iipj, 1 < j < N and the nonintersection condition of 
(£(71), . . . , £(7jv)) in the sense of Fomin (J7J, Eq.( l43|) gives the probability density function 
for joint distributions of the first passage points (x m + i9 < f l \ . . . ,x m + i$jy ) at which jj's 
first pass from the m-th chamber R{ Xm _ ltXm } to the (m + l)-th chamber R{ Xm ,x m +i} on the 
line Re z = x m , 1 < m < M. 

By the asymptotics (157)) and (155)) . the limit L — >■ 00 is taken as 

p N ( Xl ,e^;...;x M ,e W \ip) 

= lim ^(n,^;...;^,^^) 

= c N {x M ff > f%{x u o®\<p) 11 f N (x m ,e^; Xm+u e^), (44) 

given the initial condition ip G W^. 
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■ ■ ■ 















FIG. 5: Eq. ()43p gives the probability density function for joint distributions of M sets of first 
passage points on the lines Hez = x m , 1 < m < M, of iV-tuples of Brownian paths underlying in 
the nonintersecting LEBPs. 



IV. DETERMINANTAL CORRELATION FUNCTIONS 



A. Special initial condition 



Let Ax(ip) be the numerator of s\(<p) given by ([23]) . By the expansion formula fl24"]) . we 



find 



M*>) = i<( iet <N [sin^{2^ +JV - fe (cos^)^ +JV - fe + 0((cos^ +N ~ k ~ 2 )}'_ 



A' 



2^k=i( x k+ N - k ) TT s j n det 

1 1 l<j,k<N 



^COS (fj 



) x * +N - k + 0((cos< Pj ) Xk+N - k - 2 ) 



(45) 



where A\((p) is a symmetric function of {cos ¥>j}f = i with degree ^2 k=1 Afc. If we write 



7T 



/2, . . . , tt/2) as tt/2, lim^T^ i A (y>) = 5 Aj0 . Then ([22]) and ([25]) give 

f&{L,p\<p) 



lim 



2^ 2 



<^7T/2 h N {tp) K N C N (L) 

with the factor (140]) . 

Applying f[46]) to f ]39]) . we obtain the limit distribution 



h N (p) 



(46) 



^ /2 (a;,0) = lim p N {x,0\ip) = ^(h N {0)f 



)N 2 



(47) 



The fact that it is well normalized, i.e. J n/2 p^^dO = 1, is directly confirmed by using 
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the 7 = 1 case of 'the Tchebichev version' of the Selberg integral 



/:■/: 



II (Srn-tt 



l<e<m<N 



2 7 N 
3=1 



2 ,N(N-i )+ 2N TT r(i+ 7 + j 7 )(r(7j + 3/2)) 2 
+7)r( 7 (iv + j-i) + 3)' 



TT — — 
J-J- rfi 



3=0 



which is given as Eq.(17.6.4) in lOj. Under this special initial condition 7r/2, the joint 
distribution function (1441) for L — > oo becomes 



^ /2 (x 1 ,0( 1 );x 2 ,^;...;XM,0 (M) ; 

= lim p N (x 1 ,OW;x 2 ,0®;...;x M ,OM\ <p ) 



2 N Cn(xm, 



M-l 



■h N (0^) l[ f N (x m ,6^;x m+1 ,e^)h N (0 



(48) 



m=l 



for any M G N, < x\ < ■ ■ ■ < xm < oo. As a matter of course, if we set M = 1, (| un) is 
reduced to be ( 1471) . 



B. Multiple correlation functions 

For (m) G Wft, iV' m G {1, 2, . . . , N}, 1 < m < M, we put 0^ = (6[ m) , 0$), 1 < 
m < M. For a sequence {iV m }^f =1 of positive integers less than or equal to N, we define the 
(Ni, . . . , A^M)-niultiple correlation function by 

M AT j,. M 

= n n / ^" ,, pf(x 1 ,e< i ';...;x M ,e ( "»)n 7F ^ F1 T- w 

m=lj=AWi m=l 1 mh 

For 0<a;<oo,0<6><7r, we introduce two systems of functions as 

2 sm(n9) 



7r sinh(nx) 



(fi n (x,9) = 

<p n (x,6) = \ — sinh(nx) sin(n0), ra G N. (50) 

V 7T 

It is easy to confirm the equalities 

/ 4) n {x,e)H Rxl {x + i6,x' + ie')de = <f) n {x',6'), O<0'<7r, 
y H Rx ,(x + i6,x' + i9')^ n (x',9')d9' = 4> n (x,9), 0<9<n, (51) 

18 



< x < x' < oo,n 6 N. By using them, the probability density function of joint distribu- 
tions fj48l) is rewritten as follows, 



M-l 



det [^-(x!,^)] TT det [H Rx (x m + i9^\x m+1 + i9 { ™ +1) )} 

l<j,k<N K J-- 1 - l<a,/3<N x ™+i y P /J 



m=l 



x det iMx M ,e( M) )]. 



(52) 



l<p,q<N 

We have found that this product form of determinants is exactly the same as Eq.(4.5) in 
8| given for the multitime distribution function of the noncolliding Brownian motion. Then 
following the argument given in Section 4 of 8|, the following result is obtained. 



Theorem 1. Any multiple correlation function (14*91) is given by a determinant 



'■■I'M -^X., 



det 

l<j<N m ,l<k<N n ,l<m,n<M 



K'/ 2 f T f) im) -T f) {nh 



with the correlation kernel 



N 



K; /2 (x, 9; x', 9') = 0)4> n (x', 9') 



n=l 



N 



2 sinh(na;') 



n smh(nx) 

n=l v ' 



sm(nO) sinfnfl'l, if x < x' 



v 



K^ /2 (x, 9; x', 9') = Mx> WnW, d ') ~ H R* ( x> + ^ x + ^ 



n=l 



7T 



sinh(nx') 
^— ' sinh(ra) 

n=N+l V ' 



sin(n60 sm(n9'), if x > x' 



< x, x' < oo, < 9, 9' < it. 



(53) 



(54) 



(55) 



By using the terminology of probability theory, we can say that the ensemble of first 
passage points {x m + iflH : 1 < j < N, 1 < m < M} is a determinantal point process (or 



a Fermion point process) 20M22| . Note that the correlation kernel given by (|54p and (|55p 
is asymmetric in the ordering of x and x'\ K^ /2 (x, 9; x', 9') ^ ~KjJ 2 (x', 9'; x, 9) for x ^ x'. 



Such kind of asymmetric correlation kernel was first derivec 



matrix models in random matrix theory [10|, |23j . See 
Eynard-Mehta type determinantal correlations. 
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by 



30 



Eynard and Mehta for two- 



31| for recent study on the 
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FIG. 6: Conformal transformation by f(z) = e z from the domain R = {z G C : Hez > 0,0 < 
Imz < it} to the domain = {z E C : \z\ > 1, < arg(z) < ir}. The iV-tuples of Brownian paths 
in R all starting from the point m/2 are conformally transformed into the paths in £1 all starting 
from the point i. 

V. CONFORMAL TRANSFORMATION TO OTHER DOMAIN 

In the previous section, we considered a nonintersecting system of LEBPs in the half- 
infinite-strip domain, which is divided into M rectangular chambers and one half-infinite 
strip by M straight lines on Re z = x m , 1 < x\ < ■ • ■ < Xm < oo. For the underlying system 
of iV-tuples of Brownian paths, whose loop-erased parts give the nonintersecting LEBPs, 
Theorem 1 gives the determinantal correlation functions for first passage points on the lines 
Kez = x m , 1 < m < M. In order to demonstrate that the result can be conformally mapped 
to other domain consisting of a sequence of chambers in different shapes, here we show a 
conformal transformation by an entire function w — f(z) = e z . 

By this conformal transformation, the half-infinite-strip domain R = {z G R : Re 2; > 
0, < Im z < it} is mapped to the domain Q = {z = re l6 £ C : r > 1,0 < 8 < 71}. The 
rectangular chambers R{x m -\,x m }i 1 — 171 — are mapped to the chambers ^{ rm _i,r m } = 
{z = re 10 G Q : r m _\ < r < r m } with r m = e Xm ,l < m < M, and the boundary lines 
{z <E R :Rez = x m } are to the arcs of semicircles {z G Q : \z\ = r m }, 1 < m < M. 

By this conformal transformation, the paths of complex Brownian motions in R all start- 
ing from the point i7r/2 are mapped to those in Q all starting from the point i as shown by 
FigU 

The conformal invariance of the probability law of complex Brownian motions (Q implies 
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the following equality between the multiple correlation functions p^j 2 defined on R and p l N 
defined on Q, 



Al 



M 



m=l 



m=l 



where w v = ( ir 



, 4£) with w\ m) = f(x m + ]0™) = e^ m) 



{m) - '»» (m) w™) with w) m> = f(x m + i9 
1 < m < M. Since < x\ < • ■ • < xm < oo are fixed 



(m) 



g (m) 



(56) 



dw 



(m) 



diUj 



(m) 



dd 



(to) 



r m d8 



(m) 



1 < m < Af , 1 < j < N„ 



we have the following determinantal correlations for first passage points on the semicircles 
in Q. 



Corollary 2. Any multiple correlation function in Q is given by a determinant 



Mine^}^; {r^f'ft) = det 

3 3 l<j<N m ,l<k<N„,l<m,n<M 



with the correlation kernel 



1 < r,r' < oo,0 < 8,9' < n. 



N 



r ') n _ (r')~ n 



n=l 



n=N+l 



ijnTh ty — 71 



sin(n6) sin(n8'), if r<r', 
sm(n8) sm(n8'), if r>r', 



(57) 



(58) 



From (|58|) . we find that if we set r = r' > 1, 



N 



KV(re ie ,re ie '; 



— sin(n#) sin(n0') 



n=l 



_ sin((iV + 1)8) sin iV0' - sin N8 sin((iV + 1)0') 
7rr(cos# — cos^') 

for 8 7^ 8'. From it the density function p l N {re ie ) = lim £ ^ 'K 1 N (re ie , re^ 9+£ ^) is given as 

1 



(59) 



PN^ e ) 



N sin 8 - cos 8 cos(Nd) sin(N8) + sin 8 sin 2 (N8) 



iir sin 8 

For N = 3, FigJT] shows the dependence of p^re 10 ) on x = Re (re 10 ) = rcosd and y 
Im (re ) = r sin#. There are N = 3 ridges in the plots. 



(60) 
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FIG. 7: (Color online) The density function p^(re ld ) for N = 3, where x = Re (re lS ) = rcos9, 
y = Im (re 10 ) = rs'm6. There are three ridges. 

On an arc of semicircle \z\ = r > 1, < arg(z) < tt, the two-point correlation function is 
given by 

p\(re'\ re w ') = f? N (re iS )f? N (re ie ') - (K^re", re ie ')) 2 (61) 

with (J59]) and ()60]) for < 0, 0' < tt. In FigsHandEl we set r = 4 and 9 = tv/2 and plot fUD 
as a function of 9' for iV = 5 and N = 20, respectively. Due to the nonintersection condition 
for loop-erased parts, the two-point correlation function becomes zero as 9' — > 9 = tt/2. 

In general, for 1 < r < r', < 9, 9' < tt, Corollary 2 gives the two-point correlation 
function as 

4 (r'\ n — (r'\~ n 

? N (rj°, rV) = f? N (re")f? N (r'j°') + — £ 1 \ n - \i sin(n0) sin(n0') 

n=l 

x E fr A m _ (r A- m sin ( mg/ )sin(mg) (62) 

with fl60|) . For iV = 3 we set re l61 = 2e 17r ^ 2 = 2i and show by Fig{T0]the dependence of the two- 
point correlation function (1621) on x' = Re (r'e 10 ') = r' cos9' and y' = Re (r'e ie ) = r'sin^'. 
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0.5 1.0 1.5 



2.0 



2.5 3.0 



FIG. 8: For N = 5, the two-point correlation function (|6ip on an arc of semicircle \z\ = r = 4 with 
9 = ir/2 is shown as a function of 9'. There are N — 1 = 4 peaks in the plot. 




FIG. 9: For N = 20, the two-point correlation function (|6ip on an arc of semicircle \z\ = r = 4 
with 9 = 7r/2 is shown function of 9'. 



VI. CONCLUDING REMARKS 



Here we discuss the infinite number of paths limit, iV — > oo. When we take this limit in 
(!60|) . we have 

i 1 

(63) 



lim -J N (re id ) 



1 

7rr 
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y 4 



FIG. 10: (Color online) Two-point correlation function (|62p with N = 3,r = 2,9 = tt/2 is shown 
as a function of x' = Re (r'e lS ') = r' cos 9' and y' = Re (r'e l8 ') = r' sin 9'. 

That is, the distribution of the first passage point becomes uniform on an arc of semicircle. 
If we set 

r = N + u, r' = N + u', 
then the correlation kernel converges to the following as iV — > oo, 

^ N ,(N + u')e w/N )) 

2 



R- a 
9 ~N> 



n' - a 
6 ~N> 



K^it, a; u', a' 
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-(u~u')s 



sin(as) sm(as')ds, 



if u < u\ 



71 



-(u-u')s 



sin(as) sm(as')ds, if u > u' . 



(64) 



In the present paper, we have imposed special initial conditions such that all Brownian 
paths start from a single point i7r/2 for the domain R and from i for the domain Q. Study 
for general initial condition will be reported elsewhere in the future. 

At the end of this paper, we note the fact that the scaling limit of LERW is described 
by the SLE(2) path, a random continuous simple curve gen erated by the Schramm-Loewner 



evolution with a special value of parameter k = 2 29|, 
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. Kozdron 



281 ] showed that 



2x2 Fomin's determinant representing the event C(ji) fl 72 = for two Brownian paths 

le(2) ri7 = 0, where 7sle(2) and 7 denote the 



7s 



(71 ; 72) is proportional to the probability that 
SLE(2) path and a Brownian path (see also 27]]). On the other hand, Lawler and Werner 



gave a method to correctly add Brownian loops to an SLE(2) path to obtain a Brownian 
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path 341 ] . Interpretation of the results reported in the present paper in terms of 'mutually 
avoiding SLE paths' will be an interesting future problem. 
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Appendix A: Derivation of Hr l and Hqh l 



For z = x + iy e Rl, L > 0, we solve the Laplace equation 



d 2 d 2 

+ —)H Rl (x + iy, L + ip) = (Al) 



dx 2 dy 2 t 

for < p < it, by the method of separation of variables. We set Hr l (x + iy,L + ip) = 
X(x)Y(y), where description of dependence on L and p is omitted. Then we have a pair of 
ordinary differential equations 

X"(x) = cX(x), (A2) 
Y"(y) = -cY(y) (A3) 

with a constant c, which does not depend on x nor y. With the boundary condition Y(0) = 
F(tt) = 0, Eq.(jA3]) is solved as 

Y(y) = asm(ny), \fc — n, n £ N, < y < n (A4) 

with a constant a. Then Eq. (1A2j) becomes X"(x) = n 2 X(x), which is solved under the 
condition X(0) = as 

X(x) = 6sinh(nx), < x < L. (A5) 

Then we have the form 

oo 

Hr l (x + iy, L + ip) = c n (L, p) sinh(nx) sia(ny) (A6) 

n=l 
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with a series of coefficients {c n (L,p)}, where dependence on L and p is now revealed. In 
this case the boundary condition for the Poisson kernel (ITT]) becomes 

lim H Rl (x + iy,L + ip) = 5(y - p), ( A7) 

which uniquely determines the coefficients as 

. r . 2 sm(np) , , , 

c n (L,p) = - , 1 7 A8 
7r sinn (nL) 

since the Fourier series of the Dirac delta function is known as 

2 ^ 

S(y — p) — — / sin(ny) sin(np) 

71=1 

for y, p > 0. 

Following ( fl2l . the boundary Poisson kernel is obtain by taking the limit as 

H 9Rl (iy, L + ip) = lim -H Rl (iy + e,L + ip) 

2 <^ nsm(ny) sin(np) 
7r ' sinh(nL) 

n=l v ' 
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